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SUMMARY 

A probabilistic study of turbopump blades has been in progress at NASA 
Lewis Research Center for over the last 3 years. The objective of this study 
is to identify the influencing variables and quantify the effects of uncertain- 
ties in these variables on the structural response of critical structural 
components. This paper reports the effects of uncertainties in two of the 
identified variables which are geometry and material properties. The quanti- 
fied effects of the variables will provide a more reliable ground to evaluate 
the tolerance limits on the design. A methodology based on probabilistic 
approach has been developed to quantify the effects of the random uncertain- 
ties. The results of this study indicate that the uncertainties only in geome- 
try have significant effects on the structural response. 


INTRODUCTION 

Geometry and material properties are generally considered deterministic 
when designing or analyzing structural components. Contrarily, in most cases, 
these have random variations known as uncertainties in geometry and material 
properties, and sometimes even in loading. These uncertainties often conserva- 
tively are accounted for in the design by a factor of safety. It is true that 
these uncertainties are usually small in magnitude, but they sometimes have a 
significant effect on the structural performance of critical components such 
as turbine blades. Therefore, for these types of components, it is very impor- 
tant to quantify their influence on the structural response. 

The first step of this study was to develop a methodology to quantify the 
influence of these random uncertainties on the performance of structural compo- 
nents. The second step was to apply this methodology to quantify the effects 
of uncertainties in geometry, material properties, and loadings. The compo- 
nent selected for the study was the Space Shuttle Main Engine (SSME) turbopump 
blade (fig. 1). As a first phase of the study, variations associated with 
geometry and material properties have been investigated. Quantifying the 
influence of uncertainties in the loadings is underway and will be reported in 
the near future. 
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Many studies have demonstrated that probabilistic analysis methods for 
components under random loading are more reliable than deterministic 
approaches. Probabilistic methods have been predominantly used in fatigue, 
fracture mechanics, and structural reliability analyses under random vibra- 
tions (refs. 1 to 5). Particularly in fatique, improvements in estimating 
fatigue life ranged up to several hundred percents in comparison to a widely 
used deterministic approach (ref. 6). Because of the potential, the applica- 
tion of probabilistic methods have been accepted in the other fields of engi- 
neering in recent years. These fields include input loading, finite elements, 
and metallurgy. A probabilistic approach to develop a composite loading for 
SSME components is currently underway at the Battel le laboratories (ref. 7), 
and has been discussed by Shinozuka and Lin (refs. 8 and 9). Other studies 
using probabilistic finite element methods have been briefly reviewed by Das- 
Gupta (ref. 10). A variational approach for developing the probabilistic 
finite elements is currently under development by Belytschko (ref. 11). The 
usefulness and importance of the probabilistic approach, especially for turbo- 
pump blades, has been summarized by Chamis (ref. 12). 

The influence of the uncertainties in geometry and material properties in 
this study were quantified by computer simulations. The structural response 
was estimated under different levels of uncertainties using finite element 
methods. No real experiments were performed. Conducting experiments for 
studying the influence of three or more random variables with four or five 
levels and the Monte Carlo simulation technique is extremely complex and expen- 
sive, and in some cases, impossible with the state-of-the-art. Every time a 
variable or a level of variable is added to the study the number of experi- 
ments increase exponentially. 


ANALYSIS PROCEDURE 

The influence of random variations in geometry and material properties 
have been studied by modifying and using a NASA Lewis finite element code 
known as STAEBL (Structural Tailoring of Engine Blades). STAEBL modeled the 
aerofoil of the turbopump blade with 80 triangular shell elements consisting 
of 55 nodes. The blade dimensions are 1.25 by 1.00 by 0.25 in. A typical 
blade geometry is shown in figure 1. 

The influences of random variations both in geometry and material proper- 
ties, were quantified in terms of variations in the structural response of the 
blade. The structural response of the blade was characterized in terms of 
three variables. These were the first three natural frequencies, the normal 
stress (root stress) at the fixed end (junction of the aerofoil and the plat- 
form) and lateral blade tip displacements. The first three natural frequen- 
cies were considered sufficient for the purpose of methodology development. 
However, the number of frequencies can be increased to any reasonable number 
by making a simple modification in the computer code. 

Geometry uncertainties were simulated by randomly perturbing the nodal 
coordinates of the finite element mesh. All three geometrical coordinates (x, 
y, z) of each node were perturbed using random numbers generated from a com- 
puter. The generated random numbers were added vectorial ly to the nodal coor- 
dinates after the numbers were multiplied by a constant. Different constants 
were used for perturbations in different directions. The constants were 
obtained such that the maximum magnitude of perturbation in any direction 
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did not exceed 10 percent of the distance between two consecutive nodes in 
that direction. A sample of perturbed geometry is shown in figure 2. 

Material property uncertainties were simulated by perturbing the compo- 
nents of material property matrix for each finite element. Since the blade 
was very thin and the property matrix is symmetrical, only six components, 
which correspond to in-plane load resistance needed to be perturbed. Other 
components of matrix were either zero or due to symmetry did not require per- 
turbation. The maximum perturbation was once again limited to 10 percent of 
the original magnitudes of the properties. Ten percent maximum perturbation 
was justified by previous experience and a simplified spot experiment. 

To determine the extent of variables of interactions, geometry and mate- 
rial properties were perturbed simultaneously. Insignificant variables, which 
otherwise would have unnecessarily made the study more complex, were elimi- 
nated by this procedure. 

A factorial design , rather than a parametric approach, was selected to 
perform this study. This design is considered an organized and economical 
method to conduct experiments or simulations to study the influence of large 
number of variables. A large number of unnecessary tests or simulations can 
be eliminated by this design. This design helps identify insignificant varia- 
bles, thereby reducing the computation cost and complexity of the study. 
Another advantage is that this design provides the option of widening the mag- 
nitude range of the variables by allowing the additional tests outside the 
preselected ranges. The design with added tests is called Central Composite 
Design. 

A factorial design and a Central Composite Design for three variables are 
illustrated in figures 3 and 4. Solid circles on the corners, centers, and 
variable axes of the boxes in these figures represent test or simulation condi 
tions. The limitation of the design is that the results, or the models deve- 
loped using these results, are only applicable within the range covered by the 
design. A detail discussion of this method is given in reference 13. In the 
parametric approach, one variable is varied at a time. The number of tests or 
simulations needed to perform such studies using parametric approach is 
enormous . 

Simulated runs were made at each corner point of the n dimensional box 
where n stands for number of variables. For study on geometrical uncer- 
tainty, n equaled 3, for material property study, n equaled 6, and for inter 
action study, n equaled 6. The order of the runs was randomized. Also, for 
each part of study the range of each variable was extended by using central 
composite design. 


MODEL AND PROBABILITY DISTRIBUTION 

Probabilistic models were developed for all response variables for each 
part of the study using regression analysis. Coefficients and their standard 
deviations for each of the study variables were estimated. Using standard 
deviations and t-tests, it was determined whether a given variable was signifi 
cant or not. After all the models were developed, F-tests and (Chi) 
square tests were also used. These tests provided "goodness-of-f i t" estimates 
of the models and a measure of randomness of the residuals. Models which sat- 
isfied this criteria were accepted. These models are of the form; 


3 



Y = <x 0 + “i Vi + ai^2 + a 3 v 3 + + «n v n 

where Y is a response variable, o 0 , , « n are the coefficients and V] , 

, V n are the variables of the study. The residuals of the model were 

plotted and their auto correlations were plotted before a model was finally 
selected or rejected. A sample plot of residuals for one of the models is 
shown in figure 5. This plot indicates that the residuals are randomly distri- 
buted. Proabability distributions for all response variables were developed 
for selected simulated runs by replicating them a large number of times. Every 
replication employed random numbers to establish perturbation values. 


RESULTS AND ANALYSIS 

A methodology to perform probabilistic structural analysis has been deve- 
loped. It can be used to quantify the influence of random variations or uncer- 
tainties in the structural variables on the performance of a given component. 
The methodology was developed using computer simulations. These simulations 
could be used to reduce the high costs and complexity of extensive and unneces- 
sary experimentation. These advantages are further underscored by the reduced 
time an engineer/scientist must expend to arrive at these results. 

The results reported in this study, as mentioned before, are all computer- 
simulated results. The magnitudes of coefficients a], «2» - - a n for geo- 
metric variables, material property variables, and cross terms (interaction 
terms) are listed in tables I, II, and III, respectively. The significance 
test (t-test) of the coefficients in table I indicated that only variance of 
the geometry perturbations along the thickness is significant. The perturba- 
tions along other directions are insignificant and should be considered 
insignificant while setting the tolerances. The probable reason that these 
perturbations are insignificant is that they are relatively small in magni- 
tude. Although F-tests showed that all the models are good fits, they are not 
the desired models since they include insignificant variables. A model using 
the significant variable alone was also developed (not shown here). This model 
predicted estimates close to the actual values and was found to be a good fit. 

Probabilistic distributions of all the response variables were obtained 
based on 120 to 500 simulated replications. Figure 6 shows the probability 
distributions of first natural frequency as a result of geometry perturbations 
and similarly figure 7 shows the probability distributions of root stress. 

These distributions provide estimates of two items: (1) a range of variation 

of a response variable for any given random variation in the geometric varia- 
bles, and (2) probability of occurrence of a response variable. This informa- 
tion is valuable for the designer in order to adjust the design tolerance 
limits. If the range does not include any critical value, the tolerances can 
be relaxed. An obvious advantage will be that the manufacturing cost will be 
reduced. These estimates can be further improved if the number of replica- 
tions is increased. 

Significance tests were performed on the coefficients of material prop- 
erty model (table II). These tests indicated that all the coefficients are 
statistically insignificant. This means that the uncertainties in material 
properties within 10 percent of their original magnitudes did not have any sig- 
nificant influence on the response variables. The variations of the response 
variables were truly random. Therefore, controlling the material property 
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uncertainties would not necessarily improve the response. The distribution of 
all response variables corresponding to the variation in material properties 
were obtained. A sample plot of the variation in first natural frequency is 
shown in figure 8. This distribution provides estimates of range and varia- 
tions in the first natural frequency. 

Significance tests were also performed on the coefficients (table III) of 
models developed to estimate the effects of interactions. These tests indi- 
cated that all the interaction effects are statistically insignificant. This 
means that response to simultaneous variation in geometry and in material prop- 
erties are statistically independent. 

The statistical analysis was used to eliminate any insignificant varia- 
bles from the study. Since this study is continuing, these interim conclu- 
sions were used to reduce computer costs. Work is continuing in developing 
higher order models. These models include estimating variances of the response 
variables. Some interesting results have been obtained. These will be pub- 
lished in the near future. 


FUTURE WORK 

The future work includes determining the influence of random variations 
in temperature, pressure, and loadings. The loadings include random pulse 
loading, thermal loading, and random loading in the elastic and inelastic 
ranges . 


CONCLUSIONS 

1. A methodology has been developed to quantify the influence of random 
variations or uncertainties in the geometry, material properties, and loadings 
on the response variables. This methodology provides answers which may be 
used as a guideline to set the tolerance limits for design purposes. 

2. Based on this study of the turbopump blade, it was found that random 
variations or uncertainties in geometry have statistically significant influ- 
ence on the response variable. Therefore, tolerance limits should give due 
consideration to geometric variations along thickness. 

3. Random variations or uncertainties in material properties have statis- 
tically insignificant effects on the response variables. Therefore, tolerance 
on material property variations need not be precise. 

4. The methodology presented in this paper can be used for identifying 
the insignificant variables. 
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TABLE I. - PROBABILISTIC MODELS - GEOMETRIC PERTURBATIONS 


[Model: Dependent variable = Constant + Coeff * ji] + Coeff * ^ 
+ Coeff * no + Coeff * 04 + Coeff * oc + Coeff * og. F tests 
indicated that all models are good fits.] 


Dependent 
vari abl e 

Constant 

Coefficients of 

Pi 

M 2 

^3 

04 

°5 

°6 

First frequency 

6105 

-832 

-2915 

-687 

-8362 

-1897 

-9348 

Second frequency 

9475.1 

-1374 

-2961 

-408 

-11 398 

-3869 

-9095 

Third frequency 

15 792 

-9434 

27 592 

-7944 

-17 663 

-4899 

-44 147 

Root stress 

63 323 

49 707 

103 320 

4177 

88 960 

-48 191 

283 960 

Ti p di spl acement 

0.00196 


0.1287 


0.0132 


0.0067 


TABLE II. - PROBABILISTIC MODELS - MATERIAL PROPERTIES PERTURBATIONS 
(DIFFERENT BLADE PROPERTIES) 

[Model: Dependent variable = ? = Constant + Coeff ] * + Coeff£ * C ]2 

+ Coeff 3 * C 1 o + Coeff 4 * C 22 + Coeffc * C 23 + Coeff 5 ^ C 33 . F tests 
indicated that the models are not good fits.] 


Dependent 
vari abl e 
? 

Constant 

Coefficient of 

C 11 

C 12 

C 13 

c 22 

c 23 

c 33 

First frequency 

4769.7 

-239.1 

-64.1 

-55.1 

-69 

63.2 

47.1 

Second frequency 

5501 .9 

-342.6 

-25.5 

-96.7 

-250.7 

110.4 

52.2 

Third frequency 

9653 

-788.6 

- 0.8 

-61.5 

-909.6 

185.7 

-63.6 

Root stress 

97 673 

498 

-4001 

-1336 

-9962 

-12 154 

16 701 

Tip displacement 1 

0.0121 

0.0031 

0.0004 

0.0001 

0.0011 

0.0005 

0.00133 

Tip displacement 2 

0.0113 

0.0029 

-0.0006 

0.0008 

0.0008 

0.0004 

0.00133 

Tip displacement 3 

0.0012 

0.0023 

-0.0004 

0.0001 

0.0003 

0.0007 

0.0010 


TABLE III. - PROBABILISTIC MODELS - COMBINED GEOMETRIC AND 
MATERIAL PERTURBATIONS TOGETHER 

[Model: Dependent variable ? = Constant + Coeff] * O] + + Coeffg * 033 

+ Coeff 2 ] * 03 ] * C 33 . Coefficients of all cross terms (not shown) 

and material property were insignificant. Tip displacements coefficients 
to be divided by 1000 .] 


Dependent 
vari abl e 
? 

Constant 

Coefficients of 

°1 

°2 

°3 

C 11 

c 22 

c 33 

First frequency 

6185.7 

-336.2 

320.6 

-1588.1 

-60.8 

-15.37 

-44.9 

Second frequency 

8901.5 

146 

-299 

-524 

-40.1 

-68.4 


Third frequency 

15 267 

17.4 

-2307 

12 324 

-79.7 

-294 


Root stress 

85 450 

-25 410 

-39 280 

-80 

-2080 

-3260 


Tip displacement 1 

-1 .7 

- 0.10 

3.3 

-3 

- 0.1 

0.1 

- 0.1 

Tip displacement 2 

2.0 

1.5 

2.0 

1 .2 

- 0.1 

- 0.1 

0.0 

Tip di spl acement 3 

-3.1 

-3 

- 2.6 

0.70 

- 0.2 

0.1 

- 0.2 
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FIGURE 2. - PERTURBED BLADE GEOMETRY. 
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VARIABLE 3 

FIGURE 3.- FACTORIAL DESIGN FOR THREE VARIABLES. 



FIGURE 4. - CENTRAL COMPOSITE DESIGN FOR THREE VARIABLES. 
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